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We describe a minimal extension of the standard model by three right-handed neutrinos, a scalar 
doublet, and a scalar singlet (the “I'DFSZ”) which serves as an existence proof that weakly coupled 
high-scale physics can naturally explain phenomenological shortcomings of the SM. The i/DFSZ 
can explain neutrino masses, baryogenesis, the strong CP problem, and dark matter, and remains 
calculably natural despite a hierarchy of scales up to ~ 10^^ GeV. It predicts a SM-like Higgs boson, 
(maximally) TeV-scale scalar states, intermediate-scale hierarchical leptogenesis (10® GeV < Mn < 

10^ GeV), and axionic dark matter. 


PACS numbers: 12.60.Fr, 14.60.St, 14.80.Va. 

I. Introduction: The standard model (SM) and the 
paradigm of spontaneous electroweak symmetry break¬ 
ing, realised by a scalar potential 

VsM = ;r2$t$ + A($^$)2, (1) 

where « —(88 GeV)^, has proven extremely suc¬ 

cessful in explaining low energy phenomena. Neverthe¬ 
less, it fails to explain neutrino masses, the baryon asym¬ 
metry of the Universe (BAU), the smallness of the neu¬ 
tron electric dipole moment (the strong CP problem), 
dark matter, and gravity. Whether nature realises these 
phenomena in a “natural” way, i.e., in such a way that 
is (sufficiently) insensitive to physically meaningful 
quantum corrections, remains an open question. Still, 
motivated by aesthetics, the pursuit of a natural “theory 
of everything” has motivated much of modern particle 
physics. 

In the same vein, this paper describes an extension of 
the SM by three right-handed neutrinos, a scalar doublet, 
and a scalar singlet. The model can be thought of as 
an extension of the Dine-Fischler-Srednicki-Zhitnitsky 
(DFSZ) invisible axion model [HIl] by right-handed neu¬ 
trinos, and is thus henceforth referred to as the z/DFSZ. 
Notably, as we will describe in this paper, the i/DFSZ can 
explain neutrino masses, the BAU, the strong CP prob¬ 
lem, and dark matter in a calculably natural way, even 
despite a hierarchy of scales up to ~ 10^^ GeV. This 
is achieved by a seesaw mechanism, hierarchical lepto¬ 
genesis, the Peccei-Quinn (PQ) mechanism, an invisible 
axion, and a technically natural decoupling limit, respec¬ 
tively. 

The paper is organised as follows. We first detail 
the i^DFSZ, its vacuum, and its scalar sector (and con¬ 
straints). We then describe how it provides explana¬ 
tions for the strong CP problem, dark matter, neutrino 
masses, and the BAU. Penultimately, we describe our 
naturalness philosophy, identify the symmetries which 
protect each scale from quantum corrections, and study 
an example point in the parameter space. Finally we 
conclude. 


II. The vDFSZ Lagrangian: The scalar content of 
the model is a complex singlet S and two complex dou¬ 
blets $ 1,2 of hypercharge -1-1. The potential i^ 

Kdfsz = + MI2 + M^s S*S 

+ y + y + y {S*Sf 

-bAa ($ 2 ^ 2 ) +A 4 («'2®i) 

+ Ais {S*S) + X 2 S ($^$ 2 ) iS*S) 

- e$l$2S'^ - (2) 

where Mgg ^ — (lO^^ GeV)^ = —Mpg sets the PQ 
scale. Additional terms otherwise allowed by gauge sym¬ 
metry are forbidden by a global U{1)pq symmetry to be 
defined in Sec. [Vj which is essential in solving the strong 
CP problem. The e terms are necessarj0to assign a PQ 
charge to S and help to generate neutrino masses once S 
obtains a vacuum expectation value (vev). 

The only addition to the SM fermionic content is three 
right-handed neutrinos. The strong CP solution dictates 
that $1 ($ 2 ) couple to up {dp), and our solution for 
natural neutrino masses and leptogenesis requires that 
$2 couple to vp. The Yukawa Lagrangian is therefore 

-Cy = + Vu'^^iup + ydqL^2dp 

+ yeh^jep + yjL^2Vp 

+ ^yN{’^R)^Snp + H.C., ( 3 ) 

where family indices are implied, = iT 2 ^*, and J = 
2 (1) is a Type II (Flipped) j^-two-Higgs-doublet model 


^ As far as we are aware, the ^DFSZ was first discussed in Ref. [3]. 
See Refs. [US] and references therein for other minimal ap¬ 
proaches to connecting the PQ mechanism with neutrino masses. 
^ Another option is to have terms —+ H.c.) [?]• 
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(i/ 2 HDM) arrangement. We will work in the basis where 
Hn is diagonal and real. Again, additional terms are 
forbidden by the 17 ( 1 )pq symmetryj^ 

We note here that each of e —)■ 0 , i/at —)■ 0 , and 
> 0 is a technically natural limit, since they lead 
to an extra C/(l) symmetry which can be identified with 
lepton number. As well, there are two apparent techni¬ 
cally natural decoupling limits associated with enhanced 
Poincare symmetries [ 5 ]: e, Ais, A2S, j/tv —>■ 0 decouples 
S, and e, Ais, A2S, 2/1/ 0 decouples the (i>ii,S) subsys¬ 

tem. These limits will prove important in protecting the 
hierarchy of scales in the model. 

III. Vacuum: The scalar fields acquire vevs (S) = 
^'s/^/ 2 , ($i) = ( 0 ,Vi/V 2 )^. If vs > Vi, then 


Vs 



10 “ GeV. 


( 4 ) 


Expanding around this vev, the right-handed neutrinos 
acquire Majorana masses Mat = yN{S) and the scalar 
potential becomes 

K 2 HDM « mil -f TO 22 $24*2 

- + ... , ( 5 ) 

where = Mf^ + and ml2 = 

We will adopt the natural explanation of neutrino 
masses and baryogenesis detailed in Ref. [ 5 ]. This re¬ 
quires V2 ~ 0 ( 1 - 10 ) GeV achieved with mh < 0 , 771^2 > 
0 , and mf2/m22 ^ iQ Anticipating m|2 >> ^^i(A3 -I- A4), 
A2 w|, the vevs can be written 


V 2 = 


Vl 

tan P 


n2, 


m: 




22 



^ -mil V 


tan^ P 


( 6 ) 


where \Jv\Vv\ = v ~ 246 GeV and we have implicitly 
defined tan p. There is an important consistency condi¬ 
tion 2rrp22l P < \xvl to ensure m\i < 0 and avoid a 
fine-tuning for v (see Ref. i)- 

Typical values for the mass parameters are mh « 
— (88 GeV)^, 77122 10 ^ GeV, and 77122/tan ^/3 \mli\. 

Therefore, for no fine-tuning between and we al¬ 
ready expect Ais < 10 “^®, A2S < 10 “^®, and e ^ 10 “^®. 

IV. Scalar sector: The scalar mass eigenstates are, 
up to 771/77122 and 771^2/^22 corrections (see Ref. [ 3 ] for 
expressions), a CP even state [h) with m\ « Aiu^, three 
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FIG. 1. Constraints on 77122 , as labelled. Solid black contours 
are mii/GeV^ = —80^,—70^, and so on. The H/A —7 rr 
bound is taken from the CMS search m- The naturalness 
bound is only illustrative (see Ref. [^). 


heavy scalar states {H, A, H^) with masses « 77122, a PQ- 
scale neutral scalar (s) with 771^ = XsVg, and a very light 
pseudo-Goldstone boson (the invisible axion). 

Owing to the approximate t/(l) symmetry due to 
m^2/''^22 ^ 1 and v'^frrjp2 ^ the state h closely 
resembles the SM Higgsj^ In Fig. we illustrate the 
various constraints on ni22- These are the aforemen¬ 
tioned consistency condition, measurements of R —7 Xgj 
[Ullllli H/A ^ TT LHC searches (for the Type II model) 
[nm], perturbativity up to the Planck scale (TS], and 
naturalness [ 5 ]. The naturalness bound was derived in 
Ref. [ 9 ] subject to the naturalness condition we describe 
in Sec. |IX[ and we refer the reader there for details. 

V. The strong CP problem: Gauge invariance and 
renormalisability permit the physical CP violating term, 
where C is the gluon field strength tensor 
and G its dual, to be added to the SM Lagrangian. 
Bounds on the neutron electric dipole moment constrain 
6 < 10 “^® [ 12 ]. The strong CP problem is: why is 

The PQ solution [ 52 | to the strong CP problem is to 
demand a global chiral 17(1)pq symmetry; if the sum of 
the Ufj. and dp PQ charges is non-zero then, after PQ 
symmetry breaking, 0 becomes dynamical, and the vac¬ 
uum potential “selects” ( 9 ) = 0 . The pDFSZ model is 
one manifestation of this solution. 

Let us now identify the global U{1)pq symmetry. 
Defining the ?7(1)pq charge names as in Table|^ we can 
(without loss of generality) set = 0 and X^+X^ = 1 . 
Equations ([^ and ([^ set an additional six constraints on 


® In this model the right-handed neutrinos gain mass from the vev 
of S, but an alternative scenario with explicit Majorana masses 
is also possible. 

Note that, like e ^ 0 in the original Lagrangian, ^^ 12/^22 ^ 

is a technically natural limit associated with U{1)l cni. 


^ Up to v/vs corrections, at low scale the model is just the ^'2HDM 
discussed in Ref. with a very weakly coupled axion. 

^ For reviews of the strong CP problem and axions see e.g. 
Refs. 
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Field 

U{1)PQ 

Value in [Type II, Flipped] 

S 

Xs 

1 

2 

4-1 

Vi 

cos® P 

$2 

V 2 

- sin® P 


A, 

0 

Ur 

Xu 

cos® P 

Ar 

Ad 

sin® P 

II 

Ai 

f-cos®/? 

vr 

A. 

1 

4 

gr 

Ae 

!1 - 2 cos®d 


TABLE I. Charges of fields under the PQ symmetry. 


the charges, which brings the total to eight for nine un¬ 
known charges. They are completely defined by setting 
Xi = aX 2 , as long as a 1 ; it is convenient to choose 
a = — cot^ P so that the PQ current does not couple 
to the field eaten by the Z boson. The resulting charge 
values are given in Table |lj 

A final comment before moving on. In the SM, if 9 is 
set to zero by hand at some high scale, renormalisation 
implies 6 < 10 “^^ [531 [311, below the experimental 
bound. In this sense, in the SM, 0 « 0 is already natural. 
Yet this explanation remains unsatisfying, since the limit 
9 —)■ 0 cannot be identified with a symmetry. The :^DFSZ 
solution requires Ais ^ 1 , and thus one could similarly 
ask: why are the Ais « 0? At least, here, this limit is 
identified with an increased Poincare symmetry. As well, 
in the presence of CP violating new physics (such as the 
right-handed neutrinos), this solution guarantees 0 « 0 . 

VI. Dark matter: The izDFSZ axion gains a mass 


m 


a 


6 /reV 


GeV ^ 


(7) 


due to the chiral anomaly nanni, where fa ~ {S) is 
the axion decay constant, and inherits v/fa suppressed 
couplings to nucleons, photons, and electrons (for expres¬ 
sions see e.g. Ref. [T^. Stellar energy loss constrains 
/a > 4 X 10® GeV [^f] 

The axion provides a cold dark matter candidate via 
the misalignment mechanism |28H30) . wherein a signifi¬ 
cant amount of energy density is stored in coherent os¬ 
cillations of the axion field, [13] 


riah^ 


0.7 


fa 

1012 GeV 



( 8 ) 


where — tt < 0 < tt is the misalignment angle. The re¬ 
quirement that this quantity not exceed the measured 


cold dark matter energy density ~ 0.12 [3T] im¬ 

plies 


fa<2x 10“ GeV (^) ' , (9) 

with equality reproducing the observed density. If the 
PQ symmetry is broken after inflation, then the misalign¬ 
ment angle is the average value taken over many distinct 
patches, 0^ « 7 r 2 / 3 , and one obtains /a ^ 6 x 10“ GeV 
[17] Future projections of the ADMX and CAPP reso¬ 
nant microwave cavity experiments promise to probe this 
interesting region of parameter space [331IM] • 

VII. Neutrino masses: The neutrino mass matrix 
is given by 

> ( 10 ) 

where the bracketed quantity is the typical Type I seesaw 
formula |35H38j . The mass matrix is diagonalised by a 
unitary matrix U, = diag(mi, m 2 , m 3 ) = P^n, 

where m^ are the neutrino masses. It will be useful to 
express y^, in the Casas-Ibarra [33] form, 

Uy, = —VIRVI, (11) 

V2 

where P is a (possibly complex) orthogonal matrix. 

VIII. The BA U: The BAU is produced analogously 

to standard hierarchical thermal leptogenesis [40] , via the 
out-of-equilibrium, CP violating decays of the lightest 
right-handed neutrino: Vi —>■ [Sj- 

If only decays and inverse decays are considered, the 
evolution of the asymmetry is characterised (in the one- 
flavour approximation) by the decay parameter K, 


K = 


Tp 

i7|T=Mi 


mi 

m* ’ 


( 12 ) 


where F^i is the Ni decay rate, H is the Hubble rate, fhi 
is the effective neutrino mass, and m* is the equilibrium 
neutrino mass. 


Tp = 


- _ {.ylyu)iivl 
- 2M, ’ 

(13) 

Hk. 

17T2 

... 10”^ eV 

(14) 

Mpi ’ 

~ , 2a’ 

tan P 


The salient AL = 1 scatterings are electroweak, and 
those involving b quarks and (in Type II) r leptons. Since 
those rates scale with the decays and inverse decays, pro¬ 
portional to {ylyu)iii they can only represent a minor 
correction to the standard case with electroweak and t 


^ In a Type II j^DFSZ, red giants and white dwarfs constrain 

/a ^ 8 X 10® sin^ /J GeV (the white dwarf cooling fit is actually ® If inflation occurred after PQ symmetry breaking then a smaller 
improved for fa~lx 10® sin® /3) 118112511271 . 8 can be anthropically chosen, allowing fa > 10®® GeV 1321 . 
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FIG. 2. Constraints on the right-handed nentrino masses. 
The natnralness bonnd on Mn corresponds to the rough 
bound Eq. (22 1 evaluated at 17122 = 1 TeV. 


quark scatterings. The AL = 2 scatterings can how¬ 
ever constitute a significant departure from the standard 
case, particularly in FC ^ 1 scenarios dependent on ini¬ 
tial conditions. For the parameter space of interest to 
us, the generated asymmetry is safe from strong AL = 2 
scattering washout, as shown in Fig. m- As well, since 
leptogenesis in this model will be occurring at temper¬ 
atures below 10® GeV, the Yukawa couplings will be in 
equilibrium and flavour effects cannot be ignored (see e.g. 
Refs. 

These departures from the standard scenario deserve 
further detailed study. Still, we do not expect the picture 
to be dramatically changed. In particular we expect the 
Davidson-Ibarra bound [lillT] for successful hierarchi¬ 
cal thermal leptogenesis, scaled for the differing vev in 
Eq. (10), to approximately hold: 


Mni ^ 


5 X 10® GeV 
tan® P 


(15) 


This bound is depicted in Fig. 

Let us briefly demonstrate that this picture is consis¬ 
tent. A simple example configuration which achieves 
maximal CP violation and saturates the Davidson- 
Ibarra bound is (assuming normal ordering) mi <C mi 
and 


R = 


Vi - Rh 0 ^31 

0 1 0 

i?31 0 -x/l-i?ii 


(16) 


where R 31 = Kexp(i 7 r/ 4 ). Here k is related to the decay 
parameter by 


K 


0.15/A /0.05 eV^ 
tan/3 \ m 3 ) 


(17) 


and is typically 10 “® in the parameter range of in¬ 
terest. In the limit mi = 0, this corresponds to a Uyy 
which has one zero row, but is otherwise approximately 
diagonal. We note that, in this configuration, mi and 
the CP asymmetry are sufficiently stable under radiative 
corrections. The reader may check this claim against the 
renormalisation group equation (RGE) in Appendix [A} 

IX. Our naturalness philosophy: A naturalness 
problem arises when a low mass scale is subject to large 
and physically meaningful quantum corrections from a 
high mass scale. In particular, the electroweak scale can 
receive such corrections from high-scale new physics. 

The RGE formalism is a sensible way to quantify a 
naturalness problem in any perturbative quantum field 
theory. For example, in the SM, the pP{p.R) RGE is 
dominated by the top quark Yukawa, 


dpP 1 

dlnpii (47r)2 


(18) 


where yn is a representative energy scale. From the 
RGE perspective, p? is not subject to any large phys¬ 
ical corrections and can therefore be considered natural, 
in the sense that once it is set to be electroweak scale it 
stays as such within the region of validity of the model 
follows intuitively that mz) 

is not finely tuned against y?(rnz)- In the following 
we demonstrate that there exists a region of i^DFSZ pa¬ 
rameter space where our phenomenological goals can be 
achieved and the heavy PQ scale induces no naturalness 
problem, i.e., all scales remain stable under RG evolu¬ 
tion. 


X. Naturalness in the vDFSZ: Defining V = 
( 47 r)®^^jj(^ and keeping only yt,b,T,i^ Yukawas, the one- 
loop RGEs for the [Type II, Flipped] ^2HDM mass pa¬ 
rameters can be written m 


^ An additional consideration is NiNi aa annihilations. Ref. [7| 
estimates FjViAfi-^aa ~ {S) at T = which im¬ 
plies the out-of-equilibrium condition $ 10^ GeV, easily 

satisfied the parameter space of interest. 

It still may be the case that quantum gravity effects at the Planck 
scale induce a naturalness problem. However, this cannot be rig¬ 


orously computed in the absence of a cogent theory of quantum 
gravity, so an agnostic stance on this possible problem seems rea¬ 
sonable to us. In other words, we remain agnostic to the physics 
which sets the boundary condition 'f^z) $ TeV^ at 

some high scale; this is a separate problem — a hierarchy-type 
problem — which may or may not have a natural solution (see 
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Vm 

Vm 

Vm 


2 

12 


2 

11 


2 

22 


= m 


= m 


= m 


12 ~ 2 ^^ ”*" + ^Vt + 32/fe + 2 /t + Tr {ylVv^ , 

11 ~ 2^^ ”*” "*” ”*" ”*" ”*" ”*" (^•^15' + 4As) + Mgg2Xis, 

22 2 ^^ ~ 2 ^^ ”*" "*" ”*" ~*~ {yly’')^ + ™ii (4-^3 + 2 A 4 ) 

4Tr {yi^Mlf-yl) + {S)‘^X 2 S (4A2S + 4As) + Mgg2X2S, 


(19) 

( 20 ) 


( 21 ) 


where = yl^yN{S)'^ is the (diagonalised) right- 

handed neutrino mass matrix. The (S)'^ and Mgg terms 
correspond to the contribution from the heavy real sin¬ 
glet s in the broken phase. We provide the full list of 
RGBs in Appendix [X} 

These RGBs make manifest the decoupling limit 
e, Ais, A 2 S, Tr(?/j 2 /i,j/]Yj/jv) —>■ 0 which protects the scales 
from large corrections. Birst, corrections to 771^2 are 
proportional to 777 ^ 2 , reflecting the fact that e —)■ 0 re¬ 
instates a G(1 )l symmetry. Second, because the pa¬ 
rameters A 34 are reintroduced by gauge loops, mii can 
only be protected from 77722 by having 777 I 2 not too 
much larger than 777 ^]^; it was argued in Ref. [S] that 
77722 ^ few X 10^ GeV can accommodate a 10% fine- 
tuning measured at Mpi. Third, 77722 i® protected from 
the Mn by (roughly) Tr(7/fy,,7/]Y7/jv)/(47r2) < 
this translates to the sufficient condition 155] 


iWjv ^ 


3 X 10^ GeV / 77722 \ i 
tans /? VTeV/ 


( 22 ) 


for all the right-handed neutrinos, illustrated in Fig. 
Last, the are protected from the PQ scale by (again 
roughly) Xis < m‘f^/{Sy. We note that there is a lep¬ 
ton box induced correction to A 2 S; this correction is also 
bounded by m'f^/{Sy through Fq. (22). 


XI. Explicit example: As a final demonstration we 
thought it illustrative to solve the coupled set of RGBs 
for an explicit example. We consider tan /3 = 30 and 
neglect running in the following six quantities: 


Mis = - (10“ GeV)^ Mm^ = 6 x lO'^ GeV, 

{Sf = - Mls/Xs, Mn, = Mn, = 10® GeV, 

As = 0.26, yN = Mn/{S). (23) 


We have taken at the Davidson-Ibarra bound and 
below the rough naturalness bound Bq. ( [^ . We 
take T/jy according to Bqs. ([TT]), (|I^ and ([l7|) with K = 1, 


and neglect running here as well. A glance at the RGBs 
in Appendix will convince the reader that neglecting 
running in these parameters is a good approximation. 
Decoupling of the heavy degrees of freedom at mg, Mat., 
and 7772 2 should be accompanied by an associated shift in 
the Xj parameters, matching to the effective theory below 
each threshold. However, in practice, since the Xis,yu 
are so small and 77722 is not too much larger than mz^ it 
makes little numerical difference to implement this shift. 
Therefore we evolve the following parameters under the 

77 DFSZ RGBsE3 


A3(77722) — A4(77722) — 0.02, 
Ai(777z) = X2{mz) = 0.26 
yt{mz) = 0.96/sin/3, 
yb{mz) = 0.017/ cos ,5, 
yrilTlz) = 0.010/cos/3, 


Ais("7s) = 10 "^®, 

A2s(777s) = 10“^®, 

gl{mz) = 0.13, 
glimz) = 0.43, 
glimz) = 1-48. 


Their evolution is shown in Appendix 

To evolve the mass parameters we set 
77744 ( 77722 ) = —(88 GeV)^ and consider 77722 ( 77722 ) = 
0.6, 0.8,1.0, 2.0 TeV. The Ni and s are decoupled with 
step functions at their thresholds. Their RG evolution 
is shown in Fig. it is plain that the mass parameters 
in this (viable) example remain relatively small up to 
high scales, and are therefore natural according to our 
philosophy. 


XII. Conclusion : We have described an extension 
of the SM (the ‘VDFSZ”) by three right-handed neutri¬ 
nos, a complex scalar doublet, and a complex scalar sin¬ 
glet. The 77 DFSZ serves as an existence proof that weakly 
coupled high-scale physics can explain phenomenologi¬ 
cal shortcomings of the SM without introducing a nat¬ 
uralness problem. The model explains neutrino masses, 
the BAU, the strong CP problem, and dark matter, via 
a seesaw mechanism, hierarchical leptogenesis, the PQ 


e.g. Refs. l48H5ni for relevant discussion). The naturalness ques¬ 
tion induced by the high PQ scale can, by contrast, be fully 
analysed within quantum field theory, and we limit our scope to 
that. 


77 For definiteness we take a Type II arrangement, but the Flipped 
arrangement gives very similar results. 

73 We note that the parameter Ai tends to run negative, threatening 
the stability of the electroweak vacuum; nevertheless we expect 
the problem to be no worse than in the SM, i.e., we expect a 
metastable vacuum. 
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FIG. 3. Example RG evolution (see text) of 
(black solid) and (blue dashed) for 77122 = 

0.6,0.8,1.0, 2.0 TeV bottom to top. 


mechanism, and a DFSZ invisible axion, respectively. It 
contains four scales: |toii| « 88 GeV, TO 22 ~ 10^ GeV, 
Mjv ~ 10^-10^ GeV, and Mpq ^ 10^^ GeV, each pro¬ 
tected from quantum corrections by a technically nat¬ 
ural decoupling limit. The ~TeV-scale scalars and the 
invisible axion of the model will be probed in upcoming 
experiments. 
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A. RGBs : Following is the full list of RGBs in the [Type II, Flipped] model, found using PyR@TE [^. Shown 
in blue/underlined are those parameters which, for simplicity, we did not evolve in our RGB analysis. 


VMii — Mh ^6Ai — -gl — -gl + 6y/ + [0,2?/^]^ + M22 (4A3 + 2A4) + M|g2Ais, 

'DM22 = M22 ^6A2 — — -32 + ^Vb + [22 /t) 0] + 2Tr (4A3 + 2A4) + Mgg 2 X 2 s, 

'^^ss ~ ^ss (^ 4 As + Tr + MI24A2S, 

T>{S)'^ = — Tr iy^fyNj (<S')^ [i.e. the wave function renormalisation]. 


T’5{1,2,3} — {7,-3,2.3}i 

^^1 = '^9i + 29i92 + '^92 ~ -^1 (3Si + + 12Ai + 4A3A4 + 4A3 + 2A4 + 2Xlg 

+ 12X,yf-12yf+[0AXiyl-Ayt] , 

VX2 = -5i + 29i92 + ^52 ~ "^2 (3gi + 9(/|) + I2A2 + 4A3A4 + 4A3 + 2A4 + 2A25 
+ 12Ai?// - 12y^ + [4Aiy^ - Ay^, O] + 4A2Tr (yly^) - 4Tr (yly^yly^) , 

I?A 3 = -gf — -gig 2 + ^52 ~ -^3 {^9i + ^ 92 ) + + 2 A 4 ) (Ai + A 2 ) + 4 A 3 + 2 A 4 + 2 AisA 2 s 

+ A 3 ( 6 y/ + 6yl + 2yl + 2Tr (yly^)) - l^yjyl - [0,4 {ylyy)^^ yV\ , 

27A 4 = i9i92 — A 4 {^91 T T 2 A 4 (Ai + A 2 ) + 8 A 4 A 3 + 4 A 4 

+ A 4 {6yf + &yl + 2 j /2 + 2Tr (yly^)) + 12y‘^yi + [0,4 [ylyy]^^ yl] . 

27Ag = lOAg + 2 A 5 Tr (^y\iyN^ + 4 A 45 + ^A^g — 2Tr (jj\iyNylfyN^ , 

VXis = Ais 


'DX 2 S 


+ Ais i^vl + [ 0 , 2yl] + Tr , 


+ Xis ( 4 A 3 + 2 A 4 ) 


(AI) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(AlO) 


(All) 


+ A 2 S (^Vb + [2yl,0] + 2Tr (ylyu) + Tr “ 4Tr (yly^ylfyN'^ , (A12) 

27e = e (^ 2^1 ~ 2^2 + 2 A 3 + 4 A 4 + 2As + 4Ais + 4A2S + 3?// + Sy^ + + Tr (jjlyv) + Tr [y^jsjyN^ , (A13) 

T^vt = yt - ^ff 2 - 853 + ^Vt + \vl + [ 0 , yV \) > (A14) 

T^Vb = yb - 853 + \vl + \yt + [yr> o]) > (Ais) 

'Dvr = yr (^-^51 - ^yl + \yl + [3y^3y/]^ , (Ai6) 

/39 \r3 l3 1 

'^yv = yu\-^9l-^9l + ^yl + ^^{ylyv)\+ --Diag(0,0,y^)y^,0 + ^yuylyu +i^yuy^NVN. (A17) 

t^vn = ^Tr yN + yNy^NVN + VNylvu + ylylyn- (AI8) 
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B. Explicit example RG evolution : Shown below is the RG evolution of dimensionless parameters in our 
explicit example (as a function of logj^Q pLn). 


























